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Fermi satellite has detected extra spectral components in GeV energy range in several GRBs. 
Those components have power-law shapes, which may contribute to also X-ray band. The limited 
photon statistics make it difficult to determine the origin of GeV photons, namely internal or external 
shocks. We try to explain the extra components with our numerical simulations based on internal 
dissipation picture. Our leptonic model may reproduce not only the GeV excess via SSC emission 
but also the low-energy excess via the late synchrotron emission from remnant electrons. The 
hadronic models also reproduce keV-GeV power-law components by synchrotron and SSC emissions 
from secondary electron-positron pairs. In most cases the hadronic models require a much larger 
energy of protons than gamma-rays. However, the keV-GeV flat spectra detected in GRB 090902B 
is well explained with a comparable energy in protons and gamma-rays Finally, we discuss both 
advantages and weaknesses for both the leptonic and hadronic models. To overcome difhculties 
in internal dissipation models, we propose introduction of continuous acceleration similar to the 
second-order Fermi acceleration. 
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Most of the prompt emission spectra of gamma-ray 
bursts (GRBs) can be described by the well-known 
Band function [ij ; the photon number spectrum oc 
below Ep, and oc above it. The spectral peak energy 
£p are usually seen in the MeV range. One of main 
scientific targets for the Fermi satellite was to inves- 
tigate whether the spectral shape is consistent with 
the Band function even in the GeV band. The Fermi 
detected GeV photons from several very bright bursts 
(Eiso > 10^^ erg) such as GRB 080916C 0, GRB 
090902B g, and GRB 090926A 0. In such bursts 
the onset of the GeV emission is delayed with respect 
to the MeV emission. Some of them also have an ex- 
tra spectral component above a few GeV, distinct and 
additional to the usual Band function. Interestingly, 
GRB 090902B and GRB 090510 '^,^] show a further, 
soft excess feature below ^ 20 keV, which is consistent 
with a continuation of the GeV power-law component. 
While such spectra may be explained by the early on- 
set of the afterglow 0, [1| , here we pursue possilDility 
of internal-shock origins. 



II. HADRONIC MODELS 

If the spectral excesses in GeV and keV bands have 
the same origin, such a wide photon-energy range 
may imply the cascade emission due to hadrons. If 
the GRB accelerated ultra-high-energy protons, syn- 
chrotron and inverse Compton (IC) emission from an 
electron-positron pair cascade triggered by photopion 
interactions of the protons with low-energy photons 
[9IJT1I] can reproduce power-law photon spectra as seen 
in Fermi-LAT GRBs. Through Monte Carlo simula- 
tions, Asano et al. [l^] have shown that a proton 
luminosity much larger than gamma-ray luminosity 
is required to produce the extra spectral component 
in GRB 090510 as Lp > 10^^ erg s'^ (see FigH]). 
Namely, the efficiency of photopion production is very 
low. In this case, the spectrum of the GeV component 
is very hard with photon index ~ —1.6, which re- 
quires a inverse Compton (IC) contribution from the 
secondary pairs. The prominent IC component leads 
to a weaker magnetic field. This entails a lower max- 
imum energy of protons, and hence lower photopion 
production efficiency. Therefore, the required proton 
luminosity is so large in GRB 090510. 

The assumed bulk Lorentz factor in Fig[T]is 1500, 
and the emission radius is i? = 10^^ cm. If we adopt 
a smaller value of T, the pion production efficiency 
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FIG. 1: Simulated spectra with hadronic models for GRB 090510. The assumed haction of the magnetic energy density 
to the Band component photons is 10"'^, and the required proton luminosity is 200 times the luminosity of the Band 
component. 
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FIG. 2: Simulated spectra with hadronic models for GRB 090902B. The assumed fraction of the magnetic energy density 
to the Band component photons is 1.0, and the required proton luminosity is 3 times the luminosities of the Band 
component. 



would increases as tcxp/^u- oc R^^T~^. However, we 
should note that there is a lower limit to F, which 
is required to make the source optically thin to GeV 
photons. Given the photon luminosity and spectral 
shape, this minimum Lorentz factor can be estimated 
as shown in the online supporting materials in Abdo et 
al. (2009) In order to lower F, we have to increase 
the emission radius R. The lower limit of the Lorentz 
factor oc (joes not decrease drastically (since 

/? — 3, Finin oc R~^/^). The required large luminos- 
ity is rather favorable for the GRB-UHECR scenario. 



but it imposes stringent requirements on the energy 
budget of the central engine. 



On the other hand, GRB 090902B is encouraging 
for the hadronic model because of its flat spectrum 
(photon index ~ —2) ^J^. The Band component in 
this burst is an atypically narrow energy distribution 
as shown in Fig.^ f2J which may imply the photo- 
spheric emission [ij]. The hadronic cascade emission 
can well reproduce the observed flat spectra including 
the soft excess feature below 50 keV (model param- 
eters: R = 10^" cm, F = 1300). Owing to the flat 
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spectral shape of the extra component, no IC com- 
ponent is required. We can adopt a strong magnetic 
field, which enhance the photomeson production effi- 
ciency. In this case the flux of the extra component 
is relatively low compared to the Band component, 
which also decreases the required proton luminosity. 
Therefore, the necessary nonthermal proton luminos- 
ity is then not excessive and only comparable to the 
Band component luminosity. 

III. LEPTONIC MODELS 

As Corsi et al. 115| discussed, the GeV emission 
may be due to IC emission from internal dissipation 
regions. However, it seems difficult to explain spectral 
excesses in both keV and GeV bands by IC emission. 
Recently, we carry out time-dependent simulations of 
photon emissions with leptonic models In our 

simulations, as the photon energy density increases 
with time because of synchrotron emission, the SSC 
component gradually grows and dominate the photon 
field later. This late growth of the IC component has 
been observed also in the simulations of Bosnjak et al. 
(2009) [13. The resultant lightcurves show delayed 
onset of GeV emission, but the delay timescale would 
be within the approximate timescale of the keV-MeV 
pulse width. However, the longer delay compared to 
the pulse timescale such as observed in GRB 080916C 
is not explained by this effect only. 

As shown in Fig. [3](model parameters: R = 6x 10^^ 
cm, r = 1000, B = 100 G, E.^o = 10^^ erg, £e,min = 
11.3 GeV), the model spectrum obtained from our 
simulations reproduce both the low and high energy 
excesses. When the magnetic field is weak enough, 
even at the end of the electron injection, the cooled 
electrons can be still relativistic. Such cooled electrons 
continue emitting synchrotron photons. The cooling 
due to IC gradually becomes inefficient as the seed 
photon density decreases. Such late synchrotron emis- 
sion can yield the low-energy excess, while IC emission 
makes a GeV extra component. 

Another possible model is the external IC emis- 
sions [H, [Til, which can explain the delayed onset 
of GeV emission. The spatial separation between the 
source of the external photons and the site of the in- 
ternal shock region corresponds to the delay timescale. 
We also carry out a simulation for this model assum- 
ing an external emission Lgccd — 10^^ erg with 
the Band function: Epcak — 1 MeV in the observer 
frame, a = —0.6, and /3 — —2.6. The parameters 
for the internal shock are similar to those in Fig|3] as 
R = 6x10^^ cm, r = 1000, B = 100 G, Ei^o = 3x10^^ 
erg, except for ec,min = 50.6 MeV. Our simulation re- 
produces the delayed onset of GeV emission, GeV ex- 
tra component, and softening of the Band component. 
Since the external photon field in the rest frame of the 
emission region is highly anisotropic, the marginally 



high-latitude emission contributes the most to the 
flux. Thus, the simulated GeV lightcurve shows larger 
delayed onset and longer tail than that in the usual 
leptonic models (FigH]). 

IV. IMPLICATION 

The hadronic models usually require a much larger 
energy of protons than observed gamma-rays except 
for some examples such as GRB 090902B. Some may 
consider that the leptonic SSC models seem rather 
reasonable to reproduce both the GeV and keV ex- 
cesses. However, a lower magnetic field and high 
Lorentz factor, required to make an optically thin 
source for GeV photons generated via IC emission, 
leads to a very high minimum Lorentz factor (7o,min ~ 
10^) for random motion of accelerated electrons. If 
all electrons are accelerated in internal shock regions, 
such a high value may be impossible. Thus, the frac- 
tion of accelerated electrons should be small in lep- 
tonic models to explain Fermi-LAT GRBs. 

On the contrary, the minimum Lorentz factor of 
electrons 7o,min in the external IC model should be 
small to adjust the energy of scattered photons. Given 
the total energy, such small 7c,min means large number 
of electrons. If the numbers of electrons and protons 
are the same, the proton energy becomes fairly large 
(~ 1.9 X 10^^ erg in the case of FigHJ while the lumi- 
nosity of the Band component isocd = 10^^ erg s~^). 
To enhance the emission efficiency, electron-positron 
plasma should be introduced in such models. 

Despite the ability of our straightforward simu- 
lations to reproduce various observed properties of 
the GRB prompt emission, the discrepancy with the 
low-energy index a remains unexplained. The in- 
jected electrons cool via synchrotron radiation, and 
distribute below 7c.min with a power-law index —2. 
The resultant photon index becomes a ~ —1.5, while 
the observed typical values are —1.0 or harder. The 
low-energy photon index affects the photomeson pro- 
duction efficiency. We have also tried to resolve 
this problem [l^, HH considering continuous acceler- 
ation/heating due to magnetic turbulences induced 
via various types of instabilities such as Richtmyer- 
Meshkov instability The acceleration/heating 

balances with the synchrotron cooling, so the observed 
low-energy spectral index is naturally explained. Such 
effect should be included in the time-dependent code 
to reproduce the global shape of the GRB prompt 
spectra from eV to GeV. Especially, electron injection 
due to hadronic processes and succeeding acceleration 
by turbulences may explain very high 7o,min 

and GeV 

emission. 

We plan to develop the time-dependent code shown 
here to treat hadronic processes or continuous ac- 
celeration/heating. Note that the results for the 
hadronic models shown here were calculated based on 
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FIG. 3: The SSC-model fluence obtained from our time-dependent simulation (assuming z = 1). The dashed line neglects 
absorption due to EBL. 



Normalized Flux 



1 

0.8 
0.6 
0.4 
0.2 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

_,McV 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


/~\Gcv 








f^"^* 1 I ■ ■ 1 1 1 ■ ■ 1 1 ■ 1 1 ■ 1 I ■ ■ ■ 1 I ■ ■ ■ 1 ■ ■ 






1 2 3 4 5 6 7 



9 10. 



t[s] 



FIG. 4: The lightcurves for the external IC model. In the MeV band the external photons from inside region dominate, 
while GeV and eV emissions are orginated from accelerated electrons in outside dissipation region. 



the steady state approximation. We will carry out 
simulations for various situations involving dissipa- 
tive photospheres and internal or external dissipation 
or shock regions. Moreover, the code will be useful 
for simulating emissions of other high-energy sources, 
such as active galactic nuclei, supernova remnant, and 
clusters of galaxies. 
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